*Importar dataset
import excel "INSIRA AQUI O CAMINHO DO ARQUIVO DATASET", sheet("dataset survey") firstrow

*Estatistica descritiva
summarize X1 X2 X3 X4 X5 X6 R1

*Distribuicao de frequencia VD
tab R1

*Regressao logistica binominal
logit R1 X1 X2 X3 X4 X5 X6, vce(cluster CENRIO)
*logit R1 X1 X2 X3 X4 X5 X6, level(95) noconstant	

*teste de wald
test X3 X4

*AIC e BIC (modelo completo)
estat ic

*Salvar estimativas dos parametros do modelo completo
lrtest, saving(0)

*Regressao logistica stepwise forward
stepwise, pr(0.10) pe(0.05) forward: logit  R1 X1 X2 X3 X4 X5 X6, vce(cluster CENRIO)

*Odds ratio
logit, or

*likelihood-ratio test (teste de razão de verossimilhanca)
lrtest 

*AIC e BIC (modelo final)
estat ic

*Teste de Hosmer-Lemeshow
estat gof, group(10) table

*Sensibilidade do modelo
estat class, cutoff(0.6988)

*Curva de sensibilidade
lsens

*Curva ROC
lroc

*DFBETA
predict dbeta1, dbeta
scatter dbeta1 ID, mlab(ID) 
sum dbeta1

*Linktest (erro de especificacao)
linktest

*Reestimacao de Jackknife
jackknife _b, level(95) : logit R1 X1 X2 X5 X6, level(95)
